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Abstract 

The distribution of the magnitudes of seismic events is generally assumed to be 
independent on past seismicity. However, by considering events in causal relation, 
for example mother-daughter, it seems natural to assume that the magnitude of a 
daughter event is conditionally dependent on the one of the corresponding mother 
event. In order to find experimental evidence supporting this hypothesis, we an¬ 
alyze different catalogs, both real and simulated, in two different ways. From 
each catalog, we obtain the law of triggered events’ magnitude by kernel den¬ 
sity. The results obtained show that the distribution density of triggered events’ 
magnitude varies with the magnitude of their corresponding mother events. As 
the intuition suggests, an increase of mother events’ magnitude induces an in¬ 
crease of the probability of having “high” values of triggered events’ magnitude. 
In addition, we see a statistically significant increasing linear dependence of the 
magnitude means. 


1 Introduction 

The Epidemic Type Aftershock Sequences (ETAS) is a well known model in seismology 
da [la [m E]. It is a branching process in which, at each generation, each event may 
produce its own offsprings independently of the others. The magnitudes of events 
are mutually statistically independent and are distributed according to the Gutenberg- 
Richter law (5- The distribution density function of an event’s magnitude m is given 
by p{m) = 6 In 10 ■ where [3 = 6 In 10 is a positive constant 

and the cutoff parameter tuq is known as the completeness magnitude. This parameter 
is estimated from data in such a way that magnitudes exceeding it, are statistically 
in very good agreement with the Gutenberg-Richter law. This law is assumed to be 
valid both for background events and for the triggered ones and is independent on the 
magnitude of the corresponding ancestors m Recall that the background seismicity 
is the component not triggered by precursory events and is usually connected to the 
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regional tectonic strain rate; on the other hand, the triggered seismicity is the one 
associated with stress pertnrbations dne to previous shocks HU 

In this work we will use the terms “triggering event” to indicate any mother event 
that produces its own progeny [2U- Both a triggered and a background event may 
be a triggering shock. We want to investigate the distribution of the magnitude of 
triggered events, in order to assess its variation with the magnitude of the corresponding 
triggering ones. In fact, even if the correlation between subsequent events is very 
difficult to be detected and therefore is usually assumed to be absent Bn , in some 
recent works it was found statistically different from zero m m uni El EU- We expect 
that, the distribution density of triggered events’ magnitude increases (decreases) with 
the increase (decrease) of triggering events’ magnitude. Moreover, if this is true, the 
expected value of the triggered events’ magnitude should also be increasing in the same 
way. 


2 Materials and methods 

In order to obtain experimental evidence to support the above-mentioned hypothesis, 
we perform two different kinds of analysis of three Italian seismic catalogs, a Californian 
catalog and some simulated ones. The real data catalogs used here are the following. 

1. The hrst catalog includes events occurred from 16 April 2005 till 25 January 2012 
in the whole Italy. The estimated value for the completeness magnitude is 2.5. 

2. The second catalog includes events occurred in the portion of Abruzzo region 
(Italy) corresponding to the square from latitude +41.866 to +42.866 and from 
longitude +12.8944 to +13.8944. This subregion includes L’Aquila. The temporal 
interval is the same as for catalog one. The estimated value for the complete¬ 
ness magnitude is 1.8. This catalog includes the strong shock of magnitude 6.1 
occurred in L’Aquila on April the 6*^, 2009. 

3. The third catalog differs from the previous one only for the temporal interval, 
which now goes from 16 April 2005 till 05 April 2009. The value of completeness 
magnitude is estimated equal to 1.5. This catalog doesn’t include a strong shock. 

4. The fourth catalog includes events occurred from 01 January 1984 till 31 Decem¬ 
ber 1991, in the portion of the Southern California corresponding to the square 
from latitude +33.75 to +34.75 and from longitude -117.5 to -116.5. The esti¬ 
mated completeness magnitude is 2. This catalog is a portion of the waveform 
earthquake catalog relocated by Hauksson et al. in 2011 HE]. 

In all the previous cases the maximum depth considered is set equal to 40 km. 

In order to hud empirical evidence for the hypothesis that the distribution of trig¬ 
gered events’ magnitude depends on their own triggering events’ one, it seems appro¬ 
priate not to include in the catalog events that are spatially “too” distant to each other. 
This is the reason for including in the study catalogs two, three and four. More pre¬ 
cisely, we consider the third catalog to investigate the influence of a strong shock and 
the fourth one to verify the validity of our hypothesis also for a catalog relative to a 
region far away from the others and with a different seismicity. Instead, in the hrst 
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catalog there is the presence of events that are very distant to each other, as can be seen 
from the mean valnes of the distances between the events considered by us in causal 
relation, shown below. We expect that in this case the above-mentioned dependence 
becomes far less evident or absent. 

The values of the completeness magnitude have been computed with the ZMAP 
software, using Shi and Bolt uncertainty [25]. 

The simulated catalogs have been considered, instead, to demonstrate that our 
hypothesis is true also for pure simulated catalogs, not affected by any kind of “real 
effect” that may influence the analysis. The results reported here are referred to two 
simulated catalogs obtained as follows. The hrst one is simulated with the FORTRAN 
program [etasim.f] written by Ogata nanaE], just modihed in the fact that the 
number of events is random between the two input starting and ending times. We 
have chosen here the option of simulating the magnitudes with the Gutenberg-Richter 
law, instead of taking them from a given catalog. Then, we expect that there is no 
evidence of our hypothesis of conditioning. On the other hand, we simulate the second 
synthetic catalog with a program written by us very similar to the Ogata’s one, but 
adapted to our hypothesis of conditioning. More precisely, we simulate the magnitudes 
of background events again with the Gutenberg-Richter law, while the magnitudes of 
the triggered events with a new conditional probability density function with respect 
to the triggering events’ magnitudes. Since we expect that, when the triggering events’ 
magnitude m' increases, the probability of having events with “high” magnitudes must 
increase and, at the same time, the one of events with “low” magnitudes must decrease, 
we propose the following probability density function: 


p{m\m') 




1 + Gi(l - (1 - , 


( 1 ) 


where jS and a are the parameters of the Gutenberg-Richter law and the productivity 
law p{m) = respectively, and 1 > Gi > 0. We notice that this model 

reduces to the Gutenberg-Richter one for Ci = 0. For the theoretical motivation and 
properties about this formula, see |26) . 

In both the two simulated catalogs, we consider the same minimum magnitude 
value equals to 1.5 and a null learning (precursory) period. 

For each catalog, in the two types of analysis we consider four magnitudes subin¬ 
tervals contained in the magnitude range, from the completeness magnitude to the 
maximum value. The amplitude of each subinterval is opportunely chosen, for each 
catalog, in such a way in order to have a comparable number of triggered events in all 
the subintervals considered. 

We then proceed differently for the two types of analysis. In the hrst approach, we 
group the events whose magnitude belongs to each of the above subintervals. We now 
consider the events of any given subinterval. For each of them, we look at all the magni¬ 
tudes of the events occurred within a certain time interval of amplitude S* after it. We 
than group these magnitudes corresponding to all the “starting” events. Based on this 
set of magnitudes, indicated in what follows with Gm, we estimate a probability density 
function using kernel density estimation [20]. This is a non-parametric method very 
used in statistical analysis. Precisely, we consider the magnitudes (mi, m 2 ,..., m„) 
and the frequencies / = (/i,/2, • • ■, fn) of Gm- We then consider the set m of 1000 
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magnitudes equispaced from the completeness magnitude to the maximum one and we 
compute the kernel density estimator of the empirical magnitude distribution M as 


M^{m) = 



with F = 


i=l 


K 


m — nrii 


7 


( 2 ) 


where K{-) is known as kernel and the positive parameter 7 is the bandwidth [18]. The 
above formula is obtained by adapting to our case the Nadaraya-Watson kernel for 
kernel regression |24j . 

Let’s remind that a kernel is a non-negative, real-valued function such that 


K{x)dx 


1, and K{x) = K{—x) \/x G M. 


As very often done, here we use the Gaussian kernel 

/X 1 

K{x) = —;=e 2 . 

The value of the bandwidth is chosen using the leave-one-out cross-validation method, 
opportunely implemented by us ^Tj. Precisely, we consider the value that minimizes 
the quantity Ya=i \fi- fi\, where 


fi 




with 


r = J2 K 



We observe that, differently from the formula ([2]), the value ft doesn’t contribute to fi. 

The time window 5* is chosen here in such a way that two seismic events separated 
by a time larger than S* are not in causal relation. In order to determine this value, for 
each catalog we divide the whole time interval in daily subintervals. We then count the 
number of events that occur in each subinterval. Starting from this temporal sequence, 
denoted by Xt, we then compute an estimate R{6) of the autocorrelation function at 
different integer values of the time lag 6: 

.. n—S 

W) = 7-T-7> 

{n-S)V 

where n is the dimension of the sample, 6 = 0,1, 2,... and fi = ^ and V = 

are the sample mean and variance, respectively |19) . 

Then, we model R{6) by a power law model containing two parameters. These 
parameters are estimated by least squares. Finally, we hnd the value 6* such that, for 
lag values larger than 6*, the model is less than 5 * 10“^. We notice that, in the cases 
examined, this choice produces p-values always smaller than 0 . 01 . 

Due to the strong shock on April the 6 *^, 2009, the second catalog shows a clear 
non-stationary pattern. Then, for this catalog we transform the original dataset by 
considering the well-known random time change: 

X{x)dx, (3) 
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where A(-) is the ETAS rate for seismic events. The parameters of the above formula 
have been computed with the FORTRAN program [etas.f] written by Ogata, by which 
statistical inference on the parameters of the ETAS model is performed El- By this 
transformation, the process becomes stationary. 

In the second type of analysis we proceed as follows. At hrst we apply the same 
Ogata’s program as above to estimate the parameters. For each event of the cata¬ 
log, we then hnd the mother shock that most likely triggered it, by a variation of 
the Ogata’s criterion. More precisely, we consider as mother the preceding event that 
gives the highest contribute to the ETAS rate. After that, we consider the magnitudes 
subintervals strategy as in the hrst type of analysis. For each of them, we group the 
triggered events whose triggering shock’s magnitude belongs to the considered subin¬ 
terval. As in the hrst kind of analysis, for each subinterval we estimate the probability 
density function relative to triggered events’ magnitude by using the Gaussian kernel 
density estimation method described above. The value of the bandwidth is determined 
as before. 

We do not apply this analysis to the hrst catalog. In fact, it is not meaningful to 
use the pure temporal ETAS model for such a large region like the one in this catalog. 
We do not show here the results obtained for this catalog with the hrst analysis, which 
are qualitatively similar to those in the left plot at the top of Fig. |2J In fact, also in this 
case if the plot the estimated densities of the triggered events’ corresponding to four 
magnitude subintervals opportunely chosen, we can see no apparent diherences among 
the densities. As said before, this can be explained by the fact that there are many 
pairs of events that are close to each other along time, but spatially very separated. 
The elements of these pairs are erroneously put in relation in the analysis. In this case, 
the mean distance between the events of the pairs in causal relation is 140 Km. 

We implemented both methods of analysis in the MATLAB language. 


3 Results 

We present here some results obtained by the two types of the above-mentioned analysis 
for the catalogs considered. 

In Fig. [H there are the estimated densities of the triggered events’ magnitude in the 
real catalogs, obtained by the hrst and the second types of analysis (left and right plots, 
respectively). At the top, we show the results relative to the second catalog (L’Aquila 
till 2012). In this case, the spatial extension of the region analysed is far smaller than 
the one of the whole Italian catalog. We recall that, in this case, the hrst analysis has 
been applied to the dataset transformed by the random time change (jS]), due to the non¬ 
stationary pattern of the process caused by the presence of the strong shock on April 
the 6*^, 2009. The mean of the distances between the events of the pairs is about 7 Km 
and 13 Km for the hrst and the second types of analysis, respectively. From the hrst 
analysis (left plot), we notice that the increase of the referential magnitude corresponds 
to a qualitative variation of the density in agreement with our hypothesis. In fact, there 
is the increase (decrease) of the density for high (low) values of the magnitude. The 
results for the second analysis (right plot) show the same qualitative variations. The 
learning period, chosen to estimate the parameters, ends at the time of the last event 
occurred on April the 5^^, 2009. We get {fj,, K,c,a,p) = (0.304,0.06,0.104,1.57,1.39), 
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where c and p are the parameters of the Omori-Utsu law = {t + c)“^. 

In the middle of Fig. [H there are the results concerning the third catalog. In 
this case the temporal period is shorter than for the second one and ends the day 
before the strong shock on April the 6*^, 2009. The means of the distances between 
the events of the pairs are 23 Km and 17 Km for the hrst and the second types 
of analysis, respectively. Due to the absence of a very strong shock in this cata¬ 
log, the parameters used here are the averages over all the sets of parameters ob¬ 
tained by setting the learning period to 7%, 8%, ..., 20%. We get {p, K,c,a,p) = 
(0.5893,0.0219,0.0151,1.6521,1.1186). The parameters values corresponding to the 
different precursory periods show small variations from the above means. The results 
are qualitatively very close to those of the two plots at the top. This shows that our 
hypothesis of dependence is not related to the presence of a strong shock. Further¬ 
more, recalling also that the completeness magnitude is smaller for the third catalog, 
we can conclude, according to |S], that our hypothesis is not even connected to the 
incompleteness of the catalog, as instead proposed by |2]. 

At the bottom of Fig. [H we can see the results relative to the fourth catalog, that is 
the Californian one. Here, the mean of the distances is 13 Km for both the two types 
of analysis. Both from the left and the right plots, respectively obtained with the hrst 
and the second types of analysis, we get results in agreement with the above behaviors. 
It follows that, even if we analyze the events of a region in another continent, the 
hypothesis is still supported by the results of the two types of analysis. The parameters 
are here again obtained by averaging over the sets estimated for a learning period hxed 
at 7%, 8%, ..., 20%. We get {p,K,c,a,p) = (0.3729,0.0116,0.0002,0.8579,0.8879). 
Again the values obtained for the different learning periods are close to these mean 
values. 

We are going now to analyze the results obtained for the two synthetic catalogs. 

Fig. [2] contains the estimated densities of triggered events’ magnitude, obtained by 
the hrst and the second types of analysis (left and right plots, respectively), relative to 
the simulated catalogs. At the top, we plot the results concerning the catalog simulated 
with the classical Ogata’s model. In this case, the magnitudes are computed randomly 
with the Gutenberg-Richter law. According to this law, triggered events’ magnitudes 
aren’t correlated with their respective mothers’ magnitudes. This is rehected in the 
absence of variations of the densities in the four magnitude subintervals considered, 
both in the left and the right plots (hrst and second types of analysis, respectively). 
The parameters, estimated by setting the precursory at about 10%, are {p, k, c, a,p) = 
(0.62,0.02,0.013,1.72,1.11). 

Instead, a result that strongly supports our hypothesis of correlation is the one 
shown in the two plots at the bottom of Fig. |2j Here, we show the estimated densities 
of triggered events’ magnitude for the catalog simulated with our model. That is, 
the catalog in which the magnitudes are computed with the conditional probability 
density function ([T]). Once obtained this catalog, we estimated the parameters with 
the classical Ogata’s FORTRAN program [etas.f], hxing again the learning period 
at about 10%. We get {p, K,c,a,p) = (0.58,0.022,0.017,0.83,1.12). The behavior in 
these hgures is exactly the same as the one we obtain for the real data. This strongly 
supports our hypothesis. 

Finally, Fig. [3] contains the plots of the averages of triggered events’ magnitudes 
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versus triggering events’ magnitudes for both types of analysis (respectively left and 
right plots). In the two plots at the top, we consider the catalog simulated with the 
Ogata’s model and the one simulated with our model (red and black lines, respectively). 
For each of the two catalogs considered, the four triggered magnitude averages are 
normalized by the means of the average of the four values. The results of the linear 
regression analysis and the error bars are also shown. The lengths of the latter are 
given by the normalized mean standard errors. Regarding both the left and the right 
plots, one can see that there is almost no percentage variation of the triggered events’ 
magnitude in the catalog simulated with the Ogata’s model. Instead, a clear increasing 
trend is evident for the other simulated catalog considered. 

The two plots at the bottom contains the results concerning the second and the 
third catalogs, relative to L’Aquila till 2012 (blue line) and L’Aquila till 5 April 2009 
(magenta line), respectively. The means have an increasing trend in both of these 
cases. The results for the Italian and the Californian catalogs are not shown here, but 
we can say that the Californian catalog exhibits the same behavior of the plots at the 
bottom here. On the other hand, the means obtained with the first analysis for the 
whole Italian catalog are similar to those obtained for the catalog simulated with the 
Ogata’s model, as expected for the reasons explained before. 

The results obtained are statistically signihcant, as one can see from the following 
list of correlation coefficient R and p-values p. First analysis; catalog one, R ~ 0.88 and 
p ~ 0.11; catalog two, R ~ 0.99 and p ~ 0.004; catalog three, R ~ 0.96 and p ~ 0.03; 
catalog four, R ~ 0.99 and p ~ 0.0005; catalog simulated with the Ogata’s model, 
R ~ 0.61 and p ~ 0.38; catalog simulated with our model, R ~ 0.99 and p ~ 0.002. 
Second analysis: catalog two, R ~ 0.94 and p ~ 0.05; catalog three, R ~ 0.94 and 
p ~ 0.05; catalog four, R ~ 0.95 and p ~ 0.04; catalog simulated with the Ogata’s 
model, R ~ —0.92 and p ~ 0.07; catalog simulated with our model, R ~ 0.98 and 

p ~ 0.01. 

4 Conclusion 

In order to hnd evidence to support the dependence of the triggered events’ magni¬ 
tude distribution on the triggering events’ magnitude, we have applied two types of 
analysis to some catalogs, both real and simulated. In the well-known ETAS model 
the distribution of the triggered events’ magnitude is the same as that of the trigger¬ 
ing events’ one and is independent on past seismicity. Instead, our results support 
the intuitive and more realistic hypothesis of above. We notice that the probability 
density function for triggered events’ magnitude seems to vary when triggering events’ 
magnitude increases. In particular, we observe the increase (decrease) of the density 
for “high” (“low”) values of the magnitude density. This is true for all the catalogs 
with the exception of the whole Italian catalog and the one simulated by the classical 
Ogata’s model. For the hrst catalog, the absence of the variation is due to the fact 
that it contains many pairs of events temporally close to each other, but spatially very 
separated. For the simulated catalog just mentioned, the variation is absent since it 
is obtained using the standard Gutenberg-Richter law for the magnitudes. The fact 
that the catalog simulated by our magnitude model shows a clear evidence of the vari¬ 
ations, qualitatively similar to those of the real catalogs, gives strong support to the 
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validity of our hypothesis. Regarding the triggered events’ magnitude averages, we can 
see that they have always an increasing trend, again with the exception of the whole 
Italian catalog and the one simulated by the Ogata’s model. We interpret the results 
for the remaining catalogs as done before. Concluding, the law of triggered events’ 
magnitude should be a conditional probability density function that changes in shape 
with the triggering events’ magnitude. More precisely, when the latter increases, it 
may have some relative maximum for higher values of the density. It should also have 
an increasing expected value. 
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5 Figure captions 

1. Fig. [U Kernel density estimation of triggered events’ magnitude concerning the 
first and the second types of analysis (left and right plots, respectively). At the 
top, the results concerns the second catalog (L’Aquila till 2012). The considered 
intervals in which triggering events’ magnitudes fall are the following. First 
analysis, left plot: [1.8, 2.1], [2.2, 2.6], [3.3,3.8] and [3.9,5.9]; second analysis, 
right plot: [1.8, 2.2], [2.5, 3.2], [3.6,4.6] and [4.9,5.9] (in both cases, the curves 
are red, black, blue and magenta, respectively). The S* value is equal to one 
day. The optimal bandwidth value for the Normal kernel density estimation is, 
respectively for the four intervals considered, equal to: 0.22, 0.33, 0.28, 0.19 in 
the left plot and 0.25, 0.26, 0.31, 0.27 in the right one. 

In the middle, the results concerns the third catalog (L’Aquila till 5 April 2009). 
The considered intervals in which triggering events’ magnitudes fall are the fol¬ 
lowing. First analysis, left plot: [1.5,1.6], [1.7,1.9], [2.4, 2.9] and [3.1,4.!]; second 
analysis, right plot: [1.5,1.6], [1.8,2], [2.5,3.1] and [3.1,4.1] (in both cases, the 
curves are red, black, blue and magenta, respectively). The 6* value is equal to 
four days. The optimal bandwidth value for the Normal kernel density estimation 
is, respectively for the four intervals considered, equal to: 0.11, 0.14, 0.11, 0.15 
in the left plot and 0.14, 0.16, 0.22, 0.15 in the right one. 

At the bottom, the results concerns the fourth catalog (California). The consid¬ 
ered intervals in which triggering events’ magnitudes fall are the following. First 
analysis, left plot: [2, 2.2], [2.3, 2.39], [2.8, 3.2] and [3.3, 5.6]; second analysis, right 
plot: [2,2.25], [3.2,3.5], [3.5,4] and [4.6,5.6] (in both cases, the curves are red, 
black, blue and magenta, respectively). The 6* value is equal to one day. The 
optimal bandwidth value for the Normal kernel density estimation is, respectively 
for the four intervals considered, equal to: 0.11, 0.44, 0.12, 0.25 in the left plot 
and 0.18, 0.13, 0.1, 1 in the right one. 

2. Fig. [21 Kernel density estimation of triggered events’ magnitude concerning the 
hrst and the second types of analysis (left and right plots, respectively). At 
the top, the results concerns the catalog simulated with the classical Ogata’s 
model. The considered intervals in which triggering events’ magnitudes fall are 
the following. First analysis, left plot: [1.5,1.6], [1.8,2], [3.1,3.6] and [4.13,5.13]; 
second analysis, right plot: [1.5,1.65], [2,2.4], [3.3,3.8] and [4.23,5.13] (in both 
cases, the curves are red, black, blue and magenta, respectively). The S* value is 
equal to seven days. The optimal bandwidth value for the Normal kernel density 
estimation is, respectively for the four intervals considered, equal to: 0.14, 0.09, 
0.08, 0.11 in the left plot and 0.12, 0.11, 0.1, 0.1 in the right one. 

At the bottom, the results concerns the catalog simulated with our conditional 
model. The considered intervals in which triggering events’ magnitudes fall are 
the following. First analysis, left plot: [1.5,1.55], [1.7,1.8], [2.2,2.85] and [3,4.92]; 
second analysis, right plot: [1.5,1.7], [1.8, 2.1], [2.2,2.9] and [3.2,4.92] (in both 
cases, the curves are red, black, blue and magenta, respectively). The 5* value 
is equal to one day. The optimal bandwidth value for the Normal kernel density 
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estimation is, respectively for the four intervals considered, equal to: 0.05, 0.34, 
0.07, 0.18 in the left plot and 0.08, 0.08, 0.26, 0.28 in the right one. 

3. Fig. |3l Averages of triggered events’ magnitudes. At the top, the results concerns 
the catalog simulated with Ogata’s model (red) and the one simulated with our 
model (black), concerning the hrst and the second types of analysis (left and 
right plots, respectively). Regarding the left plot, the percentage means for the 
four subintervals considered for the above two catalogs are the following. Cata¬ 
log simulated with Ogata’s model; 0.9856, 0.9932, 1.0194, 1.0018 (corresponding 
to the triggering events’ magnitude 1.5506, 1.8931, 3.3302, 4.6113, respectively); 
catalog simulated with our model; 0.8256, 0.8808, 1.0508, 1.2428 (corresponding 
to the triggering events’ magnitude 1.5268, 1.7467, 2.4491, 3.4779, respectively). 
Regarding the right plot, the normalized means for the four subintervals con¬ 
sidered for the two catalogs are the following. Catalog simulated with Ogata’s 
model: 1.0052, 1.0029, 1.0011, 0.9907 (corresponding to the triggering events’ 
magnitude 1.5742, 2.163, 3.5112, 4.6786, respectively); catalog simulated with 
our model: 0.9349, 0.9827, 1.0067, 1.0757 (corresponding to the triggering events’ 
magnitude 1.595, 1.9324, 2.4614, 3.6254, respectively). 

At the bottom, the results concerns the second and the third catalog, respec¬ 
tively relative to L’Aquila till 2012 (blue) and L’Aquila till 5 April 2009 (ma¬ 
genta), concerning the hrst and the second types of analysis (left and right plots, 
respectively). Regarding the left plot, the percentage means for the four subin¬ 
tervals considered for the above two catalogs are the following. Catalog two: 
0.7460, 0.8413, 1.0748, 1.3379 (corresponding to the triggering events’ magni¬ 
tude 1.9199, 2.3491, 3.4687, 4.3342, respectively); catalog three: 0.9466, 0.9814, 
1.0233, 1.0487 (corresponding to the triggering events’ magnitude 1.5433, 1.7796, 
2.5770, 3.3818, respectively). Regarding the right plot, the normalized means for 
the four subintervals considered for the two catalogs are the following. Catalog 
two; 0.9776, 0.9844, 0.9963, 1.0417 (corresponding to the triggering events’ mag¬ 
nitude 1.9508, 2.7277, 3.9125, 5.2333, respectively); catalog three: 0.9662, 0.9801, 
0.9959, 1.0578 (corresponding to the triggering events’ magnitude 1.5433, 1.8763, 
2.7056, 3.4083, respectively). The continuous lines correspond to the results of 
the linear regression and the semi-amplitude of the error bars are the normalized 
mean standard errors. 
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